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We report comprehensive synchrotron x-ray scattering and magnetic susceptibility studies of the 
doped spin-Peierls materials Cui-ajZn^GeOs and CuGei-^SiyOa. Temperature versus dopant con- 
centration phase diagrams are mapped out for both Zn and Si dopants. The phase diagrams of both 
Cui-^Zn^GeOa and CuGei-^SiyOs closely resemble that of Cui-^Mg^GeOs, including the obser- 
vation that the spin gap is established at a much higher temperature than the temperature at which 
the spin-Peierls dimerization attains long-range order. The spin-Peierls transitions in doped sam- 
ples exhibit unusual phase transition behavior, characterized by highly rounded phase transitions, 
Lorentzian squared lineshapes, and very long relaxation times. Phenomenological explanations for 
these observations are given by considering the effects of competing random bond interactions as 
well as random fields generated by the dopants. We have also confirmed the reentrance of the spin- 
Peierls phase when the temperature is lowered through the antiferromagnetic ordering transition. 
The low temperature re-entrance of the spin-Peierls phase has been explained speculatively using a 
local phase separation scheme between the spin-Peierls phase and the Neel phase. 

PACS numbers: 75.30.Kz, 75.10.Jm, 75.40.Cx, 75.80.+q 



I. INTRODUCTION 



Low-dimensional quantum spin systems exhibit a va- 
riety of intriguing and often counter-intuitive proper- 
ties. A prominent example of such a material is the 
spin-Peierls (SP) system, which consists of an array of 
one-dimensional (ID) antiferromagnetic spin-chains with 
S=| on a deformable three-dimensional (3D) lattice. Be- 
low the spin-Peierls transition temperature, T$p, the 
spin-chains dimerize and a gap opens in the magnetic 
excitation spectrum. The discovery of an inorganic SP 
compound CuGeC>3 made possible a systematic study of 
impurity effects on SP systems.El Despite extensive exper- 
imental and theoretical efforts devoted to the study of the 
effects of doping on the magnetic and structural phase 
transitions in CuGeC>3, controversies still remain on sev- 
eral fronts. One major unresolved issue is the tempera- 
ture versus concentration phase dagr^.jjffi|bpjLh jgdthin- 
chain and between-chain dopants Bl'Mfl'SBEHO The 
phase diagram has shown certain unusual characteristics 
as revealed by different techniques: on the one hand, the 
phase behavior is so complex that there are continuing 
debates about the physics underlying nearly every sin- 



gle part of the phase diagram. On the other hand, there 
exists surprising similarities in the phase behavior for dif- 
ferent dopants, irrespective of the type and atomic size of 
the dopant ions, the magnetic moment size, or the doping 
geometry. In some theoretical models, there exists a close 
relation between the spin-Peierls transition and the sub- 
sequent antifcrrromagnetic phase transition which occurs 
at a lower temperature in samples with a dopant con- 
centration x below a critical value x c . Our own results 
suggest that the suppression of the spin-Peierls phase 
as a result of doping has no causative correlation with 
the emergence of the low-temperature antiferromagnetic 
phase except for the fact that the dopants create un- 
paired Cu spins. The lack of direct correlation is man- 
ifested most clearly in a comparative study of magnetic 
and non-magnetic ion dopings where the high tempera- 
ture suppression of the SP phase seems to follow a quite 
universal behavior with different dopants while the low 
temperature occurrence of the antiferromagnetic phase 
reflects the different magnetic properties of the individ- 
ual magnetic ion dopants. 

Currently, for within-chain doping studies, both mag- 
netic and non-magnfttifc|in^|±|a#j^^ includ- 
ing Zn, Mg and N ! .MO BliHUy iBBEy O In early 
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work on doped CuGe03, Zn was by far the most stud- 
ied dopant. A rich amount of information has been 
accumulated with a variety of experimental techniques; 
however, it has turned out that Cui-^Zn^GeOa is far 
from being the ideal system, because of concentration 
gradient effects-, and difficulties in fixing the exact Zn 
concentration.!!!] Despite these complications, we have 
continued x-ray scattering studies of Zn-doped samples 
in the hope that by directly comparing our experimental 
results with all extant ones, we can elucidate the under- 
lying physics. CuGei_j,Sij,03, on the other hand, cur- 
rently serves as the .-sole system for studying between- 
chain doping effects. LroliB Accordingly, we have carried 
out x-ray scattering experiments on CuGei_j,Si a 03 in 
parallel with our Zn studies. Finally, Mg doping, re- 
garded as the ideal dopant for within-chain dopants, has 
been studied extensively by us and oHJrKpjlte have been 
published in several previous papers .EJOHLa In this pa- 
per, we will focus on a comprehensive interpretation of 
the general phenomenonlogy for both within-chain and 
between-chain dopings. We also will present in detail a 
phenomenological model for the overall phase behavior. 

Studies of the effects of dopants on the spin-Peierls 
transition in CuGeOa started with the magnetic sus- 
ceptibility measurements by Hase et al. right after 
the discovery of. Cu.Ge03 as the first inorganic spin- 
Peierls material .Emilia In that work, Hase et al. re- 
ported a systematic study of the dopant concentra- 
tion versus temperature phase diagram using Zn as the 
dopant. Later, despite extensive experimental and the- 
oretical efforts devoted to trying to confirm this phase 
diagp m . ■-mprp-.aad- more controversies seemed to turn 
up.lt^OI£lBElEl Scattering experiments, arguably 
the most direct experimental technique for studying 
structural phase transitions, have provided invaluable in- 
sights into the phenomenology of doped CuGe03. Neu- 
tron scattering .exneiiments have been carried out by 
different groups DUlIE] In particular, Martin et alu per- 
formed the first comprehensive neutron-scattering study 
of Cui-a-Zn^GeOs. They discovered a temperature- 
concentration phase diagram in which the spin-Peierls 
phase transition temperature (Tsp) first decreased lin- 
early with Zn concentration and then appeared to show 
a plateau-like behavior above about 2% Zn doping. The 
plateau was presumed to persist up to more than 5% 
Zn doping. Nevertheless, significant broadening of the 
SP superlattice diffraction peaks was observed in a 3.2% 
Zn-doped sample which suggests that the high doping 
spin-Peierls phase may not be a true independent phase. 
This phase diagram was clearly in contradiction with 
the one obtained by Hase et al. through susceptibility 
measurements.E3 In the original phase diagram obtained 
by Hase et al., there is a linear decrease of T$p upon 
doping up to 2% Zn, which is in agreement with Martin 
et aL's results. However, Hase et al. reported the disap- 
pearance of the SP transition above 2% Zn doping; this 
was indicated by the disappearance of the characteristic 
dip in the magnetic susceptibility curve, which reflects 



the opening of the SP spin gap and therefore ordinarily 
signifies the onset of the SP transition. 

At first glance, the above two phase diagrams seem to 
be irreconcilable with each other. However, if one takes 
a closer look at the experimental specifics, one finds that 
the discrepancy originates in part from the different ap- 
proaches of two experimental techniques. Susceptibility 
measurements, as a macroscopic experimental technique, 
are primarily sensitive to long range order and less so to 
local fluctuations. Neutron scattering, on the other hand, 
given sufficiently coarse resolution and long scan times, 
can easily reveal the existence of short range fluctua- 
tions. The neutron scattering results reported by Martin 
et al. strongly suggest that the high doping SP phase 
is of short range order and hence not easily detected by 
magnetic susceptibility measurements. Thus, care must 
be taken in determining the phase diagram of the high 
doping phase. Similar ambiguities have been discovered 
in studies of Si-doped samples. 

Thus, a systematic and detailed study using both sus- 
ceptibility and scattering techniques should help to re- 
solve some of the extant discrepancies. In this paper, we 
report a detailed synchrotron x-ray scattering and mag- 
netic susceptibility study of the spin-Peierls transition 
in both Cui_ x Zn x Ge03 and CuGei_ y Si y 03. The paper 
is organized as follows: First we present the results for 
Cui_2;Zn x Ge03 as a prototypical example to elaborate 
on the major experimental results. The general conclu- 
sions are applied to our study of CuGei_ y Si y 03. We 
then present a phenomenological interpretation of our 
combined results, in which the important roles played by 
competing interactions, random fields and phase inhomo- 
geneity will be emphasized£3 

In the next section, we describe the sample preparation 
and the experimental techniques used in this study. In 
Sec. Ill, our experimental results from x-ray scattering 
and magnetic susceptibility measurements are presented, 
including the experimentally determined phase diagrams 
for both Zn- and Si-doped CuGeOs. Finally, we will 
discuss various theoretical and experimental issues in Sec. 

ivl 



II. EXPERIMENTAL DETAILS 

A. Crystal Growth 

The single crystals used in the experiment come from 
two sources. The CuGei_j,Si a 03 single crystals were 
grown at the University of Tokyo using the traveling 
solvent floating zone method and preliminary sample 
characterization and magnetic susceptibility results have 
been published elsewhere. E£l The impurity concentration 
of these samples has been determined by inductively 
coupled plasma atomic emission spectroscopy (ICP-AES), 
with an accuracy of about 0.1%. The Cui_ x Zn x Ge03 
crystals were grown by the same method at MIT. Sam- 
ples have a typical mosaic spread of less than 0.05° full- 
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TABLE I: Characterization of Cui-zZn^GeOs crystals. Tsp 
and Tjv are determined from the magnetic susceptibility mea- 
surements (see Fig. 1). 
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width at half-maximum at major Bragg diffraction peaks. 
The typical sample size used for both x-ray scattering and 
magnetic susceptibility measurements was about 3x3x1 
mm 3 . We prepared our samples by cutting them from 
the end of regular growths which have lengths over 5 
cm. This procedure effectively eliminated the gradient 
induced in crystal growth caused by the different doping 
levels of the seed crystals. 

In order to differentiate the intrinsic physics from ex- 
traneous effects caused by simple concentration gradi- 
ents, a detailed and careful characterization of the sam- 
ples is crucial. The electron probe microscope analysis 
(EPMA) method was used to determine the Zn doping 
concentrations. Twenty evenly spaced spots covering the 
whole sample surface were used in the EPMA measure- 
ments. The variations of concentrations were recorded as 
the uncertainties in the sample concentration. The sam- 
ple characterization results are summarized in Table I. As 
shown in Table I, the actual Zn contents are always lower 
than the nominal ones in the high doping range, which is 
consistent with previously reported results.t_Mj However, 
despite the great disparities between these two concentra- 
tions, the concentration gradients and fluctuations found 
over the surfaces remain relatively small and compara- 
ble to those in the Mg-doped samplesB2l Moreover, the 
concentration variations have similar magnitudes for all 
samples, which suggests comparable effects, if any, on the 
observed SP phase transitions. 



B. Synchrotron x-ray scattering 

High resolution synchrotron x-ray scattering is a pow- 
erful tool in studying structural phase transitions. High- 
flux synchrotron x-ray radiation enables us to study both 
the critical fluctuations and the long range order. Our ex- 
periments were carried out at MIT-IBM beamline X20A 
at the National Synchrotron Light Source. The 8.5 keV 
x-ray beam was focused by a mirror, monochromatized 
by a pair of Ge (111) crystals, scattered from the sample, 
and analyzed by a Si (111) analyzer. Carefully cleaved 
samples were placed inside a Be can filled with helium 



heat-exchange gas and mounted on the cold finger of a 
closed-cycle refrigerator. The measurements were carried 
out in each sample at the (1.5, 1, 1.5) SP dimerization 
peak position with the (H K H) zone in the scattering 
plane. The high momentum resolution of the synchrotron 
x-ray beam enabled us to measure intrinsic correlation 
lengths as large as 5000 A; any system with a correlation 
length larger than this was considered to possess long- 
range order (LRO). 



C. Magnetic susceptibility 

The magnetic susceptibility of the identical samples 
used in the x-ray experiments has been measured using 
a commercial SQUID magnetometer (Quantum Design 
MPMS). The samples were mounted with the c axis par- 
allel to the applied magnetic field, and the data were 
taken in a magnetic field of 500 Oe. 



III. EXPERIMENTAL RESULTS 
A. Magnetic susceptibility — Zn doping 

Figure |l| shows the magnetic susceptibility as a func- 
tion of temperature for all of the Cui-^Zn^GeOs sam- 
ples. The overall features are similar to previously re- 
ported results, QBE3£3 showing a rapid suppression of the 
SP phase transition temperature upon doping and the 
appearance of a low temperature antiferromagnetic (AF) 
phase. For doped samples, the kink anomaly, which is a 
characteristic of the SP transition, is rounded. We de- 
termine Tsp in the susceptibility measurements from the 
maximum of the derivative of the kink anomaly and the 
Neel temperature, TV, simply from the low-temperature 
cusp temperature. For samples with £=0.0235, 0.0255, 
0.34, and 0.38, a sharp low-temperature cusp which sig- 
nifies the onset of the AF state is observed. On the 
other hand, the susceptibility curve of the x=0.027 sam- 
ple shows a broad low temperature peak with two cusp 
positions possibly corresponding to two Neel tempera- 
tures. 

This is reminiscent of similar results obtained for Mg- 
doped samples by Masuda and coworkers £3 Those au- 
thors carried out susceptibility measurements on a series 
of Cui^Mg^GeOa samples, and found a double peak 
feature in the low temperature susceptibility data, for x 
in the neighborhood of 0.023. In their study, it was ar- 
gued that there exists a first order phase transition as 
a function of x, between a dimerized antiferromagnetic 
phase and a uniform antiferromagnetic phase, so that 
around the critical concentration x c ~ 0.023, two AF 
transitions appear because of the microscopic phase sep- 
aration between these two phases. Detailed comparison 
of the magnetic susceptibilities of Zn- and Mg- doped 
samples yield many similarities, as expected. In Ref. 
O, the sample with the highest Mg-doping without the 
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FIG. 1: The magnetic susceptibility as a function of temper- 
ature for various Zn-doped samples. Each curve is shifted by 
10" 3 emu/g for clarity. 



broadening feature had a Neel transition around 2.9 K 
and XMg = 0.0229. The low threshold started with 
XMg = 0.0237 in which there existed a broadened cusp 
spanning from 3.2 K to 3.9 K. The cusp peaked more 
to the low temperature side with the weight gradually 
shifting to thciiigh temperature side as the concentration 
was increased.E3 We note in Fig. ^ similar behavior in our 
Zn-doped samples, with the x = 0.0255 Zn-doped sample 
showing a single cusp at 2.9K, while the x = 0.027 sam- 
ple shows a broadened cusp over the temperature range 
from 3.1 K to 3.8 K. For x=0.027, the left side of the 
cusp is much higher than the right side suggesting that 
this sample is closer to the lower temperature threshold of 
the peak broadening phenomenon. Unfortunately, we do 
not have a range of doping concentrations from £=0.027 
to x=0.034 which could give us a more complete view of 
the gradual shift of the weights. The important point is, 
however, that we have confirmed that the regime where 
the susceptibility shows significant broadening coincides 
with the region where, as we shall discuss in the next 
subsection, the SP phase begins to lose LRO while sub- 
stantial superlattice scattering intensity is nevertheless 
still observable Jii Therefore, the Neel transition broad- 
ening must somehow be correlated with the loss of LRO 
of the SP phase. 



B. X-ray scattering — Zn doping 

The instrumental resolution function was experimen- 
tally measured at the (3, 0, 0) Bragg peak, which has 
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FIG. 2: Longitudinal scan profile of the superlattice peak, 
(1.5, 1, 1.5), at T=5 K for a;=0.038. The solid and dashed 
line shows the correponding results fitted by Lorentzian and 
Lorentzian squared lineshapes convoluted with the instrumen- 
tal resolution function. 



a Q value close to that of the (1.5, 1, 1.5) dimerization 
peak. The determination of the intrinsic line-shape is 
complicated by the existence of the two-length scale phe- 
nomenon as has been discussed in previous work by some 
of us.E-3 Since two superimposed profiles corresponding 
to two length scales were always present in the criti- 
cal scattering in the close vicinity of the SP transition 
for x < 0.02 samples, a unique determination of the in- 
trinsic cross-section was difficult. However, only a single 
scattering profile was observable in the critical scattering 
for x > 0.02 samples, where dopant effects play a more 
important role than the near surface dislocation effects 
which are presumably responsible for the large length 
scale fluctuations in samples with LRO. For typical sec- 
ond order phase transitions, the critical scattering cross 
section assumes a Lorentzian lineshape. Any deviation 
from this line-shape can provide information on the ef- 
fects of disorder on the underlying phase transition. 

In Fig. |[ we show the low temperature (1.5, 1, 1.5) 
superlattice peak profile for the x=0.038 sample. It is 
evident that the resolution function is much narrower 
than the intrinsic cross section. Both Lorentzian and 
Lorentzian squared cross sections have been used to fit 
the data. The Lorentzian squared lineshape gives a dis- 
cernably better fit than a simple Lorentzian lineshape. 

Figure ^|(a) shows the temperature dependence of the 
(1.5, 1, 1.5) dimerization peak intensity for various 
Cui-^Zn^GeOs samples. We normalize the intensity of 
each sample data set using the dimerization peak inten- 
sity at 5 K. As is obvious in the figure, in sharp con- 
trast to the transition in pure CuGeOa, the SP transi- 
tions for the doped samples are noticeably rounded, and 
further the rounding increases progressively with increas- 
ing dopant concentration. We believe that this substan- 
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FIG. 3: The (1.5, 1, 1.5) SP peak intensity (a) and the 
corresponding inverse correlation length (b) as functions of 
temperature for various Zn-doping concentrations. The in- 
set shows the impurity-concentration dependence of the low- 
temperature SP inverse correlation length 



tial rounding is intrinsic and cannot simply be accounted 
for by local concentration gradients for the following rea- 
sons: (i) The impurity concentration-is too low to cre- 
ate any substantial correlated effect t3 (ii) There exists 
a monotonic increase of the rounding in the order pa- 
rameter as the doping concentration is increased, while 
experimentally the measured concentration gradients are 
of similar magnitude for all of the samples studied. For 
example, the a;=0.038 sample has a concentration vari- 
ance of less than 0.2%, from Table I. However, that sam- 
ple's SP dimerization peak intensity increases sufficiently 
gradually as the temperature is lowered that it does not 
achieve saturation for the whole temperature range stud- 
ied. In sharp contrast, no significant rounding was found 
in the x=0.011 sample which has a similar concentration 
variance. 

Not only does the order parameter temperature depen- 
dence change dramatically upon doping, the correspond- 
ing saturation correlation length also behaves quite dif- 
ferently from that in the undoped sample. Figure ||(b) 
shows the temperature dependence of the longitudinal in- 
verse correlation length for the different Cui-^Zn^GeOa 
samples. The minimum inverse correlation length which 
the sample achieves as a function of x is plotted in the 
inset. Unlike the pure material, in which very near T$p 
the inverse correlation length decreases linearly to zero 
as Tsp is approached from above because of the second 



length scale critical fluctuations, the inverse correlation 
length of the doped samples appears to approach zero in 
an asymptotic way without a clear signature of the exact 
temperature at which it becomes zero. This complica- 
tion makes an unequivocal determination of T$p difficult. 
Within our experimental accuracy, samples attain LRO 
(that is, the correlation length exceeds 5000 A) only for 
x < 0.027. Substantial intensity is still observable at the 
superlattice peak positions for the x=0.034 and x=0.038 
samples which have x > x c = 0.027, although the su- 
perlattice peaks are not resolution-limited in width. The 
critical dopant concentration, x c , which we have deter- 
mined for Cui-zZn^GeOa is slightly higher than that for 
the Cui-^Mg^GeOs samples. However, this could be due 
in part to the different techniques used in determining the 
concentrations . 

In Fig. [|, we show the results from the fits for the 
peak intensity, which for LRO is the order parameter 
squared, and the inverse correlation length, £ _1 for sam- 
ples with £=0.011 and x=0.0235. Also plotted is the dif- 
fuse scattering intensity at the wing (Sq = 0.022 A -1 ) of 
the dimerization superlattice peak, which is a measure of 
the amplitude of the critical fluctuations. The transition 
is slightly broadened as compared to the SP transition in 
pure CuGeOa. This is reflected both in the smearing of 
the order parameter onset temperature and the temper- 
ature range of the diffuse scattering. The inverse correla- 
tion length follows more or less a linear relationship with 
temperature. For the samples with a;=0.0235, the broad- 
ening in the temperature dependence of both the order 
parameter and the diffuse scattering has increased signifi- 
cantly. The corresponding inverse correlation length first 
exhibits a linear relationship as a function of tempera- 
ture with decreasing temperature. However, close to the 
transition, T < 9.5 K, the rate of the divergence slows 
down and £ _1 appears to approach zero asymptotically, 
that is, with zero slope. It is interesting to note that if 
the linear behavior is extrapolated to intersect the tem- 
perature axis, one obtains a temperature which roughly 
matches that of the maximum of the diffuse scattering. 

SP transitions in doped samples have shown unconven- 
tional characteristics that are absent in typical second- 
order phase transitions. One of these is that when the 
doped system undergoes a SP transition, the times for 
the system to reach phase equilibrium are anomalously 
long. In x-ray scattering experiments, this is reflected in 
a continuous change in the peak intensity and width as 
a function of time. This anomaly manifests itself in pro- 
nounced thermal hysteresis in the order parameter mea- 
surements if the system is not kept sufficiently long at 
each temperature. This slow dynamics behavior is not 
observed in pure CuGe03. Under the same experimen- 
tal conditions, upon rapid change of temperature by sev- 
eral degrees, the pure system achieves equilibrium within 
minutes, in sharp contrast to the highly doped samples. 
We characterized the slow dynamics by monitoring the 
intensity of the superlattice peak after a sudden change 
in temperature as a function of time while keeping the 
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FIG. 4: Top: Peak intensity of the superlattice peak (empty 
symbols) and the diffuse scattering intensity at (1.5 1.03 1.5), 
plotted as a function of temperature. Bottom: Inverse corre- 
lation length of the critical scattering along the a, b, and c 
directions. 



final temperature constant. The time scans for the sys- 
tem in equilibrium show an oscillating behavior with the 
centerline remaining flat. The oscillations reflects the 
temperature fluctuations at the sample place and agree 
well with the temperature fluctuation patterns recorded 
from the temperature controller. 

Figure |B| shows representative time scans for :r=0.0235 
for both heating and cooling temperature runs on a semi- 
logarithmic scale. After the sample was quenched from 
8K to 7K, it was observed that the peak intensity in- 
creased monotonically with time and did not saturate 
after more than half an hour. There are, in addition, 
oscillations identical to those observed in the pure sam- 
ple. This reveals that the slow dynamics are not caused 
by poor thermal contacts. The slow dynamics behavior 
exists even upon change of temperature by an amount as 
small as 0.05K as shown in Fig. |](b). The linear rela- 
tionship in Fig. U suggests that the peak intensities are 
following time logarithmic behavior. In addition to the 
anomalous behavior of the peak intensity, the peak width 
also changes continuously. A slow increase of intensity is 
accompanied by a slow decrease in peak width and vice 
versa. This means that the time dependence originates in 
slow growth or contraction of the domains after a sudden 
change in temperature. 

To understand the origin of the slow dynamics, we first 
recall that for SP systems, the phase transition occurs as 
the neighboring spins pair up and create a local gap in 
the excitation spectrum, and eventually LRO is achieved 
when the 3D phase coherence of these local dimeriza- 
tions is established. In pure CuGeC>3, these two events 
happen at the same temperature. For doped samples, 
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FIG. 5: (a) Semi-logarithmic plot of the peak intensity of the 
superlattice peak as function of time after a sudden temper- 
ature change from 8.0 K to 7.0 K for 1=0.0235. (b) Similar 
plot after a temperature change from 8.03 K to 8.07 K. 



on the other hand, we observed substantial superlattice 
intensity even when the system had only SRO, which sug- 
gested that SP order had been established locally though 
the system had difficulty in achieving global phase co- 
herence. Similar behavior has been, observed in random 
field Ising model (RFIM) systemsBj The slow dynamics, 
then, results from the expansion or contraction of the 
SRO region as a function of time. However, we do not 
have quantitative information on the time evolution of 
the peak width. 



C. X-ray scattering — Si doping 

As discussed in Sec. I, one of the findings in the stud- 
ies of doped CuGeOs is the close resemblance of the 
phase diagrra m for , b oth within-chain and between-chain 
dopings.l3fl't3'E3E3Ea That is, substitution on the Cu site 
and substitution on the Ge site have yielded very simi- 
lar effects on the SP phase, including rapid suppression 
of the SP transition temperature, occurrence of a low- 
temperature AF phase, and the concomitant co-existence 
of the SP and AF phases. 

In this section, we describe synchrotron x-ray diffrac- 
tion experiments on CuGei_j,Sij,03 samples. As ex- 
pected, most of our experimental observations are sim- 
ilar to the Zn-doping results presented in the previous 
section. However, we find that Si-doping is two to three 
times as effective in destroying the SP phase as compared 
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to Zn- or Mg- doping, in agreement with earlier workers.H 
A representative set of longitudinal (parallel to the 
scattering vector) x-ray scans around the (1.5, 1, 1.5) SP 
dimerization peak position for the y=0.0019, 0.0104, and 
0.0171 samples at various temperature is shown in Fig. ^. 
Similar to the Zn-doping case, the low-temperature SP 
peak width is resolution limited in the y=0.0019 sample, 
indicating that the SP phase possesses LRO. On the other 
hand, the SP peaks in the y=0.0104 and y=0.0171 sam- 
ples are broadened at all experimentally accessible tem- 
peratures, and therefore only SRO is present. To extract 
the correlation length, the data for all samples were fit 
to a two-dimensional convolution (longitudinal and trans- 
verse directions) of the experimentally measured Gaus- 
sian resolution function with the intrinsic cross section 
taken as Lorentzian squared. 

In Fig. ^|, the peak intensity and the associated inverse 
correlation length obtained from the fitting are shown 
as functions of temperature. For the y=0.0019 sample, 
with less than 0.2% of Si substitution, the SP transi- 
tion temperature has been reduced by 1 K. However, 
other than that, the £=0.0019 sample behaves analo- 
gously to pure CuGeC>3. The effects of doping can be 
more obviously seen as we gradually increase the dopant 
concentrations. For the samples with higher doping lev- 
els, smearing of the order parameter onset temperature, 
which is reminiscent of the behavior in Cui-^Zn^GeOa 
and Cui^zMgzGeOa, can be identified. In addition, over 
the whole temperature range studied, the y=0.0104 and 
the y=0.0171 samples never achieve LRO and the corre- 
lation lengths at low temperatures saturate at the finite 
values of about 300 A and 90 A, respectively. Moreover, 
for the y=0.0104 sample, this saturation of the corre- 
lation length at low temperatures is preempted by the 
low temperature re-entrance of the SP phase, which is 
characterized by a suppression of the dimerization peak 
intensity and a concomitant decrease in the correlation 
length. The re-jeatrance, as already elaborated in our Mg 
doping studies ,ta can be understood as arising from local 
phase separation of the SP and AF micro-regions. For 
the y=0.0171 sample, despite the large error bars due to 
the low intensity, the reentrance at around 4 K is still ob- 
servable from the drop in the peak intensity upon cooling 
from 4.5 K to 4.1 K. This temperature range also includes 
the Neel temperature of this sample. 
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FIG. 6: Representative longitudinal scans at the (1.5, 1, 1.5) 
SP dimerization peak position. The solid lines are the results 
of fits to a Lorentzian squared lineshape as described in the 
text. The instrumental resolution function is shown as dashed 
lines. 
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D. Phase diagram 



FIG. 7: Peak intensity (top) and inverse correlation length 
(bottom) of the SP peak as a function of temperature for the 
CuGei-^SiyOa samples. 



In order to construct comprehensive temperature ver- 
sus impurity concentration phase diagrams from our x- 
ray scattering and susceptibility measurements, we define 
two characteristic temperatures, both related to the SP 
transition but defined by the physics on different length- 
scales. The higher transition temperature T m , is defined 
as the temperature corresponding to the peak tempera- 
ture of the critical fluctuations as measured from diffuse 
scattering intensity. We defined the counterpart of T m 



in the bulk susceptibility measurements as the tempera- 
ture at which the derivative of the susceptibility (dx/dT) 
reaches a maximum. The lower transition temperature 
T s is defined as the temperature at which the lattice 
dimerization achieved LRO, that is, when the inverse 
correlation length reaches zero, which in our case means 
< 0.0002 A" 1 . 
For the x=0 data shown in Fig. pi there is no am- 
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biguity in the definition of T$p. The super lattice peak 
intensity vanishes above T=14.1 K while the correspond- 
ing correlation length diverges precipitously at the same 
temperature. In the inset, we show the diffuse scatter- 
ing intensity at the wing (Sq = 2.1 x 10 _3 A _1 ) of the 
dimerization peak. It is clear that the critical fluctuation 
amplitude reaches a maximum at 14.1 K as well. The 
susceptibility measurement shown in the bottom panel 
of Fig. U also exhibits a sharp kink anomaly at 14.1 K, 
at which temperature the derivative shows a sharp peak. 
From these observations, we find Tgp = T m = T s = 14.1 
K for x=0. Physically this means that the temperature 
at which the spin gap appears coincides with the temper- 
ature at which long-range coherence of the dimerization 
is achieved. This perfect agreement by two different tech- 
niques, one magnetic and one structural, proves that the 
SP transition is a well-defined second order phase transi- 
tion and it also demonstrates the equivalence of the two 
definitions for the transition temperature. However, as 
discussed in our previous papers and as is evident from 
the data presented in this paper, a discrepancy develops 
between these two definitions as one progresses into the 
doped regime. It is then difficult to define T$p unam- 
biguously. 

In Fig. ^, we show the experimental results for the 
x=0.0235 sample. In the upper panel where x-ray mea- 
surements of the dimerization peak are presented, we find 
smearing of the onset temperature of the order parame- 
ter, with substantial intensity starting to appear around 
10 K while the critical fluctuation amplitude reaches a 
maximum around 9.3 K. At the same time, the cor- 
relation length does not diverge until the temperature 
reaches T s 7 K on cooling. In the bottom panel, mag- 
netic susceptibility measurements are shown on the same 
temperature scale. The kink anomaly corresponding to 
the SP transition is rounded and the derivative of the 
susceptibility shows a broad peak around T m w 9.5 K. If 
a linear extrapolation method, which was adopted in Ref. 
|l2| , is used, we obtainTsp ~ 10 K. These seemingly con- 
tradictory results partially explain the difference among 
the varied phase diagrams reported in the literature. 

In Fig. |Io|, we plot the T — x phase diagram for 
Cui-zZn^GeOa using the above definitions from both 
x-ray and susceptibility results. Our phase diagram is 
in several ways distinctive: The phase boundaries de- 
termined by different definitions of T$p do not match 
each other except for x=0, with the divergence increas- 
ing as the doping level increases. The phase boundary 
determined from T s is nearly vertical around x=0.027. 
No LRO can be attained above this concentration level, 
although substantial intensity at the dimerization peak 
position is still observed. Note that the AF transition in 
the susceptibility measurement at a:=0.027 is also broad- 
ened. On the other hand, the phase boundary deter- 
mined from T m follows an approximately linear variation 
with the Zn concentration. In particular, for x=0.034, 
while T s cannot be defined due to the lack of LRO, one 
can still determine T m because there still exists substan- 



CuGeO, 




FIG. 8: Top: Temperature dependence of the peak intensity 
(empty) and the inverse correlation length (filled) of CuGeOs. 
Inset: Temperature dependence of diffuse scattering intensity 
at 5q = 2.1 x 10~ 3 A -1 . Bottom: Temperature dependence 
of the magnetic susceptibility (empty circles) and its temper- 
ature derivative (filled circles). 



tial intensity at the dimerization peak position in the 
critical temperature regime. 

In Fig. [n], we-also present the T-y phase diagram for 
CuGei-ySiflOaS The same definitions for T s and T m 
are used. The phase diagram is very similar to those 
of Cui-zMg^GeOa and Cui-^Zn^GeOa, with a linear 
suppression of the SP phase upon doping and a critical 
concentration around a; c =0.009, above which no LRO can 
be observed. The value of x c rj 0.009 is slightly less 
than half that of the critical concentration x c =0.021 and 
x c =0.027 in the Mg- and in Zn-doped cases, respectively. 
Also, one should note that the slope of the linear line 
in Fig. [ll] is two to three times that in Fig. [To], which 
suggests that the Si-doping is more than twice as effective 
in destroying the SP phase as the within-chain Mg and 
Zn dopants. Previous workers have rationalized this bv. 
noting that the Si-dopant affects more than one chain.tj 
In addition, as we shall discuss in the next section, the 
Si will create a random elastic field on the Cu site which 
in a pseudo-spin description will act like a random field 
which is also destructive of the order. 
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FIG. 9: Top: Temperature dependence of the peak inten- 
sity (empty) and the inverse correlation length (filled) of 
Cuo. 9765 Zno. 0235 Ge03. Also plotted in square symbols is the 
temperature dependence of the diffuse scattering intensity at 
Sq — 2.2 x 10 -2 A -1 . Bottom: Temperature dependence of 
magnetic susceptibility (empty) and its temperature deriva- 
tive (filled). 
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FIG. 10: Experimentally determined T — x phase diagram of 
Cui-zZn^GeOs- 
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FIG. 11: Experimentally determined T-y phase diagram of 
CuGei-ySi H 03. The solid line shown in Fig. [l(j is rescaled to 
illustrate the y = 2x and y = 3x scaling. 



IV. DISCUSSION 

The first conundrum in the diluted CuGeC>3 problem 
is why the SP phase is destroyed by only ~ 2% Zn or Mg 
and even less in the Si case. In comparison, a much higher 
doping concentration is required to destroy the Haldaiwj 
gapped state in the S — 1 spin chain.E2l Khomskii et a/.cJ 
first proposed a soliton model as a starting point for the 
study of the doped CuGeOa problem. The solitc 
is consistent wit.b-.hoth numerical simulations^!! 
and experimentsEj'Eij However, Khomskii and coworkers 
note that if diluted CuGeC>3 is treated in a simple per- 
colation picture, then the soliton model naturally yields 
a critical concentration x c ~ 10%. Here x c is determined 
by the concentration at which the average dopant spac- 
ing equals the soliton width. This implies that there must 
be a separate mechanism that is responsible for the rapid 
destruction of the SP phase. 

In order to discuss this, we follow Harris et al£3 and 
map the dimerized SP system onto an effective 3D Ising 
model, in which the two dimer configurations possible in 
a given chain are associated with the up and down states 
of the Ising pseudo-spins as shown in Fig. |l2|(a). When a 
few Cu ions are replaced by Zn ions in an isolated chain, 
it is energetically favorable for the Ising pseudo-spin to 
change sign across the Zn dopant. In Fig. |l2|(b), we show 
that the dopant effectively creates an antiferromagnetic 
bond in the ferromagnetic pseudo-Ising model. Note that 
the dopants can cut the chains into segments with cither 
an even or an odd number of spins in them. The "even" 
segments favors the dimer configuration in which all the 
spins are paired, since this configuration has lower energy 
in comparison with the other one, which has two un- 
paired spins at the ends of the segments. For the "odd" 
segment, on the other hand, the two dimer configura- 
tions are energetically equivalent. According to numer- 
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FIG. 12: (a) Pseudo-Ising model mapping for single chain 
without dopant and (b) with dopant (c) Pseudo-Ising model 
with interchain interchain 



ical simulations. El however, the soliton in the odd seg- 
ments prefers to be in the center of the segment, thereby 
creating two smaller segments with opposite pseudo-Ising 
spins. 

As temperature is lowered, a uniformly ordered phase, 
as shown in Fig. [l2|(c) is preferred, since the interchain 
interaction prefers to keep all the dimers in phase in order 
to achieve 3D LRO. Therefore, depending on the energet- 
ics of the interchain interaction and the soliton creation 
energy, locally ordered dimers with the "wrong" pseu- 
dospin may be nipped at low temperatures. We note a 
similar mapping by Mostovoy et al.fB in which the equiv- 
alent fields were taken to infinity to solve the model. 

The essential physics thus involves a competition be- 
tween the interchain interaction which favors a uniform 
phase of the dimerization and the local soliton energy 
which is minimized by a change in phase of the dimer- 
ization at the dopant site. It is this competition be- 
tween intrachain and interchain energies which leads to 
the fragility of the LRO spin-Peierls state. 

Another system of this kind is a magnetic system with 
both ferro- and antiferromanetic bonds. For this system, 
the ordering problrap-jwith competing interactions has 
been well addressed. c3o In particular, Aeppli, Maletta 
and coworkers nmpo.sed a, phenomenological cluster spin- 
glass mode]c3cd'c3'L§Oc3 to explain the evolution of the 



ground state from ferromagnetic LRO to spin-glass. 

We should emphasize that our observation in doped 
CuGe03 of a Lorentzian squared cross section instead of 
Lorentzian is significant. By way of contrast, in a sim- 
ple dilution induced percolation problem£3 even above 
the percolation threshold, the diffuse scattering should 
always be Lorentzian, which is characteristic for fluctu- 
ations of thermal or geometrical origin only. The occur- 
rence of a Lorentzian squared lineshape implies coronet- 
ing interaction or random field dominated physics.c!3 In 
the Si-doped case, because of the low symmetry of the 
Ge site the dopant will generate both competing inter- 
chain interactions and random fields. Microscopic calcu- 
lations are required to determine which are more impor- 
tant quantitatively. 

In addition to the unusual SP transition, the AF tran- 
sitions in doped CuGeOs samples also exhibit abnormal 
characteristics. Neutron scattering experiment on 3.4% 
(nominal) Zn-doped CuGeOsEJ revealed that the AF or- 
der parameter did not saturate even at 1.4 K. In ad- 
dition, the transition itself was considerably rounded. 
Fitting the order parameter using a Gaussian distribu- 
tion of the transition temperature, Hase et al. obtained 
(3 = 0.22 ± 0.02 which did not fit intp-aay extant 3D uni- 
versality class. In the /iSR studies ,Ol£i the Neel order- 
ing process was described as an inhomogeneous process 
with islands of freezing spins emerging from the param- 
agnetic phase. Although the transition was mistakenly 
interpreted in a spin glass picture, it revealed that the 
Neel transition was not a homogenous process despite 
the fact that it yields a 3D LRO magnetic state. We 
believe that the inhomogeneous ordering of the AF tran- 
sition is due to the coexistence of an inhomogeneous SP 
phase and the AF phase involving local phase separation. 

In order to understand the competition between the 
SP and AF phases, it is illuminating to compare 
Cui^Zn^GeOs with Cui-^Ni^GeOs. In the case of Ni- 
doping casej] due to the non-zero magnetic moment of 
the Ni ion, the AF ordering takes place at a quite differ- 
ent temperature from that of Zn- or Mg- doping, while 
the effect on the SP phase boundary seems to be the 
same. This experimental observation suggests that the 
SP phase transition and the AF transition should be 
treated independently, that is, the suppression of the SP 
ordering can not be simply attributed to the competition 
from the AF ordering. Although these two phases are 
intricately correlated as evidenced by the re-entrance of 
the SP phase upon the onset of the AF phase, the on- 
set of the AF ordering is only indirectly related to the 
destruction of the SP phase. 

As discussed in detail in the section on Si-doping and 
in Ref. |l8|, we observed at approximately the same tem- 
perature as the onset of the AF ordering, that the SP 
phase underwent a reentrant transition which was evi- 
denced by the increase in the dimerization peak width 
upon cooling through Neel temperature. Additionally, 
the reentrant process was accompanied by slow dynam- 
ics and hysteresis. The reentrance and the hysteresis are 
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FIG. 13: (a)The free spin configuration above the AF transi- 
tion and (b) Antiferromagnetically ordered islands below the 
AF transition. 

unusual in the sense that they coincide with the onset 
of the AF ordering over the doping range studied, which 
implies that the reentrance is directly related to the AF 
ordering. 

As shown in Fig. [H?], we propose a local phase separa- 
tion scenario to explain this unusual AF transition and 
the coexistence of the AF and SP phases. As shown in 
Fig. [l^, solitons bear double identities in the system, they 
involve both structural anti-phase domain walls for the 
SP phase and the free spins for the underlying magnetic 
system. For materials with x < x c or y < y c , the temper- 
ature is decreased below the SP transition temperature, 
the interchain interactions lead to the creation of a LRO 
SP phase. For even segments, two solitons are created 
at the end - see Fig. |L2](c), while for odd segments, the 
solitons that were located at the center will move to the 
end of the segment to bind to the dopants. As a re- 
sult, a LRO SP order in fact creates many solitons that 
bind to the dopant sites. These free spins have antifcr- 
romagnetic correlations among them, and the dimerized 
spins will be partially polarized to mediate the magnetic 
interactions. At the Neel temperature, the spins even- 
tually develop AF order. We believe that the key factor 
in inducing the reentrance is the non-uniform nature of 
the AF state. Above the AF transition, the free spins are 



disordered. However, as the system is cooled, the antifer- 
romagnetic interaction energy will become relevant. The 
free spins are actually mobile and can move freely in each 
segments. However, due to the confining potential of the 
interchain interations the solitons become bound to the 
dopants. The onset of the AF state, however, will com- 
pete with this configuration because, the lowest energy 
state will be achieved by moving all these solitons to- 
gether, in other words, phase separation into soliton-rich 
and soliton-poor regions, which will inevitably disrupt 
the originally established SP order. 

The final state is the result of this competition between 
the SP order and AF order. Since the solitons are bound 
by the dopants, the phase separation can happen only 
locally. It is not currently known whether the resulting 
state involves islands or stripes or some other geometry. 
Further theoretical work must be done in this direction. 
Clearly, however, the slow dynamics and re-entrance can 
naturally be explained in this way. 
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